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INTRODUCTION 


A side-looking imaging radar can obtain good azimuth resolution 
either by utilizing a long narrow beajn antenna, or by developing a 
"synthetic aperture" which requires a time integration of the phase of 
the scattered signal. The problem with using real apertures from space 
is that resolution of a few tens of meters requires an antenna that is 
a few hundred meters long. The problem with the synthetic apertvire is 
that the ocean waves are moving as the synthetic aperture is developed. 
The latter problem is the subject of this paper. 

The statistical theory of rough surface scattering has been used 
with success to define the average microwave scattering cross-section of 
the ocean. Such a statistical theory; however, is not appropriate for 
interpreting radar imagery since the scene is deterministic. Therefore, 
in order to analyze the synthetic aperture radar, an ocean-like surface 
must be generated, and the calculations of the magnitude and phase of the 
scattered energy must be undertaken on that particular surface. . 

In order to derive a closed form solution of an ocean wave image, it 
is assumed that the ocean wave propagates only in one direction and that 
the wave is monochromatic. The antenna is assumed to illuminate the 
stirface in the near nadir direction, and geometric optics is assumed to 
govern the scattering process. The image is reconstructed by doppler 
ordering of the phases of the scattered signal. Parametric calculations 
are presented as a function of integration time, platform altitude to 
velocity ratio, ocean wavelength and electromagnetic wavelength. It is, 


; 
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shown that unless the system parameters are judiciously selected serious 
distortion of the image can occur. 

Discussion of Figures 

Figure 1 shows the two-dimensional geometry for the imaging of 
specular points by means of doppler ordering the back- scattered electric 
field as the antenna beajn moves over the specular points. In this figure, 
R is the radius vector to a point along the X-axis, H and v are 
the altitude and velocity of the moving platform and Xj^ is the position 
of the ith stationary specular point. The time reference is chosen so 
that at t = 0 h is the perpendicular distance to any given ith specular 
point . 

The phase relationship for the return from the ith specular point 

is given by the equations in the figure, where the latter is an approxi- 
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mat ion for large h. Note that the phase contains a t term due to the 
variable doppler shift of the specular point as it passes through the 
antenna beam. 

Figure 2 shows the simplified system model used in this paper for the 
analysis of the scattered return. The received signal is mixed with 

a reference chirp to re-order the doppler components. The signal then 
goes through a low pass filter to produce S^, whose phase depends 
linearly on time. This signal is then integrated over time T and 
squared to produce the intensity, The resultant signal is a 

series of sinc^ functions of width X = h/2vT, where X is the 
resolution of the image. For analytical simplicity, a zero IF stage 
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which passes both positive and negative frequencies has been postulated. 

Figxire 3 shows a very simplistic model of a moving ocean wave that 
is being imaged ♦ The wave is a single harmonic gravity wave of wave 
number K = 2 jt/L (L = wave length) and radian frequency Cl =V^ ^ 
where g is the acceleration due to gravity. In this geometry^ the 
antenna beam 9g is assumed to be small enough so that only those 
specular points at normal incidence are illuminated. Therefore, reflec- 
tions occur at the peaks and troughs of the ocean wave; i.e., a double 
harmonic is imaged. Also the amplitude of the return from each specular 
point is the same: only an uninteresting phase shift of 90 ^ occurs 

between the peaks and troughs. 

Because of the motion of the wave, the specular points will move 

in time in proportion to the phase velocity v^^ - Vg/K , as given by 

the equation in the figure. Since X - (t) is squared additional terms 
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proportional to t will appear in the exponential. The time integration 
of signal Sj will therefore be Fresnel integrals instead of simple sinc^ 
functions. 

The details of the derivation explicitly give 


where 
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As the antenna scans off nadir^ the specular points will move 
towards the even zero crossings of the X-axis, so i;hat the illumination 
approaches a spacing equal to L, the ocean wavelength. Indeed a simple 
calculation shows that the specular points are identically spaced L 
apart when the 0, the viewing angle, is 


0 


X -1 / 2nA . 

9 = tan ( ) 

c L 


There are no specular points for incidence angles greater than that given 
above. Therefore, the specular points occur at distances ^ apart at 
nadir; and approach a spacing of L at the critical viewing angle 0^ , 
beyond which, they entirely disappear from view. There is no unique 
relationship between the amplitude and length of an ocean wave; however, 
the rms slope (A/L) generally increases with increasing wind speed 
and decreasing depth of the bottom. Under these conditions, some specular 
returns will occur with reasonable probability for viewing angles up to 
20 or 50 degrees. At greater angles of incidence, Bragg scattering will 
dominate. 

Figure 4 shows the radar image of a "frozen" ocean wave (Vpj^ = O). 

The value of 40 sec. for the h/v ratio is approximately the value for 
a high flying transonic jet aircraft. Intensity peaks occiir every 20 
meters, or one-half an ocean wavelength. This figure shows that the 
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resolution of the image improves as the integration time increases. 

Figures show how the image of a moving specular point is distorted 
as the integration time increases. The calculations assume that the ocean 
wave propagates with a phase velocity v^^^ • For. T = 0-5 sec. 

(figure 5)^ the intensity distribution closely resembles the sinc^ function. 
As T increases to 1.0 sec. (figure 6), the width of the central peeth is 
only slightly greater than the width of the sinc^ function. However, 

large side lobes are beginning to appear. When the integration time 
increases to 1.2 sec. (figure J), the pair of side lobes are almost as 
intense as the central lobej and severe image distortion occurs. Finally,- 
for T = 1.1 seconds (figure 8), destructive interference causes a minimum 
to occur at the specular point position. 

Figure 9 shows the distortion of the specular point image as a function 
of the electromagnetic wavelength A or h/v for L = lO m and T = 1.0 sec. 
With h/v = lO sec., the image degrades as the electomagnetic wavelength 
becomes shorter than 0, 5 m. A re-scaling of parameters shows that with 
A = .25 m. , the distortion is reduced as h/v increases. However, an 
increasing h/v will result in a larger value of i.e., poorer resolu- 

tion. 

In figure 10, T = 1.0 sec., A =0.25 m and h/v = 40 sec., and L is 
allowed to vary. As the ocean wave increases, the phase velocity also 
increases; and the image is distorted. On the other hand, as the ocean 
wavelength becomes smaller, the specular points are closer to each other 
and the image of a given specular point smears into the image of its 
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neighbor. This figure therefore implies that there is both a minimum 

and a maximum ocean wave that can be imaged. 

Figiires 11 and 12 outline methods by which the imaging of an ocean 

wave can be optimized. Figure 11 indicates how a distorted image may be 

reconstructed by hardware implementation. If the Chirp frequency is 

tuned by an amount equal to 1 - 2Vpjj/V; then the mixing and filtering of 

signals and results in a signal S^, which contains no quadratic 
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time varying phase terms. The integrated signal then becomes a sine 
function, and the specular point is resolved to within x^. The diffi- 
culty is that one must know what the ocean wavelength is in order to 
measure it I Another difficulty is that this method will restore only one 
wavelength of what is generally a spectral ensemble of wavelengths. 

The other method (figure 12) is to make an appropriate selection 
of the electromagnetic wavelength to image the range of ocean wavelengths 
of interest. This can be done by utilizing the two conditions given in 

the figure. The first criterion x^ < L . A resolves the minimum 

o — min' 

ocean wavelength as reviewed from near nadir. The second criterion 

> v^^ T follows from inspection of the Fresenel integral expression 

given earlier. Basically, the latter inequality states that a given specular 
point must remain Within the resolution cell during the integration period. 
Otherwise image distortion will result. These criteria are conservative 
because the specular points will tend to be spaced a full wavelength apart 
as the viewing angle increases; and the calculations show that latter 
criterion can be "slightly" violated before serious image distortion occurs. 
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If the inequalities are replaced by equalities, and if T is 
eliminated, one obtains the result in the box. If h/v is fixed, this 
equation states that the greatest range of ocean wavelengths can be 
illuminated by choosing the shortest practical electromagnetic wavelength. 

CONCLUDING REMARKS 

Geometric optics has been used as a scattering model in order to 
interpret the synthetic aperture image of a moving ocean wave. The results 
basically show that the image of a monochromatic ocean wave will be distorted 
if the product of the integration time and the phase velocity of the wave 
exceeds the resolution of the synthetic aperture. 

When the monochromatic wave is viewed in the nadir direction, specular 
points (images) occur at increments of one-half an ocean wavelength. As 
the viewing angle increases, the images occur at one-wavelength intervals 
and then disappear vjith increasing viewing angle. 

Numerical techniques can be used to generate a surface which spectrally 
resembles an ocean wave, and calculations performed within 20 to 50 ° of 
nadir should be physically meaningful. At steeper viewing angles, the 
model must include the Bragg Scattering Mechanism. 



IMAGING RADAR GEOMETRY 



IMAGING RADAR SYSTEM 



Figure 2. 









IMAGING OF OCEAN WAVES 



A COMBINATION OF FRESNEL INTEGRALS 


Figure 5 • 


IMAGE OF A STATIONARY WAVE 



Figure . ^ . 



IMAGE OF MOVING WAVE 

T = 0. 5 sec 


IMAGE OF MOVING WAVE 

T = 1 sec 


INTENSITY 




Figure 5- 


Figure 6 . 


IMAGE OF MOVING WAVE 

T= 1.2 sec 



IMAGE OF MOVING WAVE 



Figure 7 . 


Figure 8 , 


IMAGE VERSUS X OR h/v 



Figure 9 . 


IMAGE VERSUS OCEAN WAVELENGTH 



T = 1 sec 
X = 0.25 m 
h/v = 40 sec 


Figure 10. 



SYSTEM OPTIMIZATION 

A. HARDWARE IMPLEMENTATION 



Figure 11. 


SYSTEM OPTIMIZATION 

B. SELECTION OF PARAMETERS 
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Figure 12. 







